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Abstract 



^ In this paper, we propose a new scheme for anisotropic motion by mean 

-^1- curvature in R^. The scheme consists of a phase-field approximation of the 

CTn motion, where the nonlinear diffusive terms in the corresponding anisotropic 

r ■ Allen-Cahn equation are linearized in the Fourier space. In real space, this 

^^ corresponds to the convolution with a kernel of the form 

in 

. . We analyse the resulting scheme, following the work of Ishii-Pires-Souganidis 

^ on the convergence of the Bence-Merriman-Osher algorithm for isotropic 

motion by mean curvature. The main difficulty here, is that the kernel i^^^t 

is not positive and that its moments of order 2 are not in L^(R^). Still, 

C^ we can show that in one sense the scheme is consistent with the anisotropic 

mean curvature flow. 

1 Introduction and motivation 

In the last decades, a lot of attention has been devoted to the motion of interfaces, 
and particularly to motion by mean curvature. Applications concern image pro- 
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cessing (denoising, segmentation), material sciences (motion of grain boundaries 
in alloys, crystal growth), biology (modelling of vesicles and blood cells). 

1.1 Motion by isotropic mean curvature 

The simplest case of motion by isotropic mean curvature concerns the evolution 
of a set ftf C M^ with a boundary d^t of codimension 1, whose normal velocity 
Vji is proportional to its mean curvature hc 

Vn{x) = /^(x), a.e. X G F^, (1) 

with the convention that hc is negative if fi^ is a convex set. It at t = the initial 
set r^o is smooth, then the evolution is well-defined until some time T > when 
singularities may develop [2J. 

Viscosity solutions provide a more general framework, that defines evolution 
past singularities, or evolution from non-smooth initial sets. If ^ is a level set 
function of f^o, i.e., 

f^o = |x G M"^ ; g{x) < o} , dQo = {x G M"^ ; g{x) = o} , 

and if u denotes the viscosity solution to the Hamilton- Jacobi equation 

u(0,x) = g(x), 

then the generalized mean curvature flow fit starting from Qq is deflned by the 
0-level set of i/ [211 EDI m El 



nt 



{x E M"^ ; u{t, x) < O} , dn{t) = {x G M"^ ; u{t, x) = o} . 



Alternatively, one can deflne the motion by mean curvature as the limit of 
diffuse interface approximations obtained by solving the Allen- Cahn equation 

du . 1 



dt 



Au-^W\u), (2) 



e 



where e is a small parameter (that determines the width of the diffuse interface) 
and where W{s) = ^ ^"^ is a double well potential. This equation can be 
viewed as a gradient flow for the energy 

Mu)= I ('-\Vu\'' + -W{u)] dx. 



Modica and Mortola ^3 ES] have shown that J^ approximates (in the sense of 
F- convergence) the surface energy cw J where 

J(n) = / Ida and cw = [ \/2W{s)ds. 

JdVL Jo ^ 

Existence, uniqueness, and a comparison principle have been estabhshed for ^ 
(see for example chapters 14 and 15 in [2J and the references therein). 
Let Ue solve ^ with the initial condition 



where d{x, Q) denotes the signed distance of a point x to the set Q and where 
the profile q is defined by 

q = argminj^ (^l'\s) + W^(7(s))) ds;^ G i?L(M), 7(-oo) = +1, 

7(+^) = -l, 7(0) = ^}. 

Then, for smooth motion by mean curvature JUII], or for generalized motion 
by mean curvature without fattening [31 [21] , the set 

approximates Q{t) at the rate of convergence 0(e^|logep). 

The Bence-Merriman-Osher algorithm [9j is yet another approximation to 
motion by mean curvature. Given a closed set E C M^, and denoting xe its 
characteristic function, one defines 

ThE^LeR"^ ; u{x,h)>^Y 

where u solves the heat equation 

f^{x,t) = Au{x,t), t>0 X G M^ 
\u{x,0) = Xe{x). 

Setting Eh{t) = T^^/^^E, where [t/h] is the integer part of t/h, Evans [2D], and 
Barles and Georgelin [4J have shown that Eh{t) converges to Et^ the evolution by 
mean curvature from E. 



1.2 Motion by anisotropic mean curvature 

We use the framework of the Finsler geometry as described in [8j . Let : M^ ^ 
[0, +oo[ denote a strictly convex function in (7^(M^\{0})), which is 1-homogeneous 
and bounded, i.e., 

\ A|e| < m < A|CI C e K^ 

for two positive constants < A < A < +00. We assume that its dual function 

(/)° :R^ ^ [0, +00 [, defined by 

</>''(r) = sup{r.e; </'(e)<i} 

is also in C^(R^\{0})). Given a smooth set E and a smooth function ix : R^ ^ R 
such that dE = | x G M^ ; u{x) = [, we define 

• the Cahn-Hoffman vector field n^^ = (/)?(Vi^). 

• the (/)-curvature hicf) = div(n(^). 

We say that E(t) is the evolution from E by ^-curvature, if at each time t, the 
normal velocity Vn is given by 

Vn = -i^cf^n^- 

As in the case of isotropic flows, one can define motion by ^-curvature using a 
level set formulation, i.e., following the level lines of the solution to the anisotropic 
Hamilton- Jacobi equation 

ut = (f)%Vu) 0^^(Vi/) : V\. (3) 

Existence, uniqueness and a comparison principle have been etablished in [TBI ESI 

The anisotropic surface energy 

J{n) = / (j)\n)da. 
Jdfi 

can be approximated by the Ginzburg-Landau-like energy 

J,,^(n) = / (U%Vuf + -W{u)) dx, 



and its gradient flow leads to the anisotropic Allen-Cahn equation [T] 

du . 1 



dt 



A^u-^W'{u). (4) 



e 



The operator A^^ := div ((/)2(Vix)0^(Vi/)] is called the anisotropic Laplacian. 

The Bence-Merriman-Osher algorithm has also been extended to anisotropic 
motion by mean curvature. One generalization was proposed by Chambolle and 
Novaga [12j as follows: Given a closed set E, let Th{E) = | x G M^ ; u{x, h) > ^\, 
where u{x,t) is the solution to 

-rrix.t) = A4,u(x,t), t>0 x E M^ 
at 

^u{x,0) = Xe{x). 

Define then Eh{t) = T^ ^E. The convergence of Eh{t) to the generalized 
anisotropic mean curvature flow from E is established in |12| . The result holds 
for very general anisotropic surface tensions and even in the cristalline case. How- 
ever, because of the strongly nonlinear character of A^^, the numerical resolution 



of (1.2) is much harder than in the isotropic case. 

Another generalization of the Bence-Merriman-Osher algorithm has been stud- 
ied by Ishii, Pires and Souganidis ^27j. The main idea is to represent the solution 



u of (1.2) as the convolution of xe with a geometric kernel. More precisely, Let 
/ : R^ — ^ R be a function which satisfles the following conditions 

(Ai) Positivity and symmetry : 

f(x) > 0, f{—x) — /(x), and / f{x)dx — 1 

(^42) Boundedness of the moments : 

/ \x\'^ f{x)dx < +OC, 

< / (1 + \x\^)f{x)dn'^-^ < 00, for all p G S^-\ 
{A3) Smoothness : 

p ^ /(x)(iH and p ^ XiXjf{x)dT-r~^ are continous on S' 
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Given E C M^, define ThE = {x G M^ ; u{x, h)>\^, wliere 

u{x,h)= / kh{y)xE{y -x)dy, 

with the kernel 

Kt{x) = ^f{Vix), xgM^. 

They showed [27] that T^^E converges to the set E{t) obtained from E as 
the generalized motion by anisotropic mean curvature via the Hamilton Jacobi 
equation 

ut = F{D^u,Du) 

where 

F{X,p) = (^1^ f{x)dn''-\x)^ (^-^ 1^ {Xx , x) f{x)dn''-\x] 

This result raises a natural question: Given an anisotropy 0^, can one find a 
kernel /, so that the generalized fromt dE{t) defined by the associated Hamilton 
Jacobi equation evolves by (/)-mean curvature ? This problem has been addressed 
by Ruuth and Merriman [33j in dimension 2: They propose a class of kernels and 
study the corresponding numerical schemes, which prove very efficient. However, 
their appraoch cannot be generalized to higher dimensions. In contrast, our 
algorithm is not specific to the dimension 2. 



1.2.1 A new algorithm for motion by anisotropic mean curvature 

In this work, our objective is to extend Ishii-Pires-Souganidis' analysis to study 
the following algorithm. Starting from a bounded closed set E C M^, we define 
an operator T^E by 

ThE^^x^R'^- u{x,h)>]^, (5) 

where u solves the following parabolic equation: 

\«(x,0) = xe{x). 
Denoting by F{u) the Fourier transform of a function u, 

T{u){i) = / u{x)e-'^^'''< dx, 
6 



the operator A^^ is defined by 

A^n = ^-i(-47rV(0'-^(^)(0)- 

and can be seen as a linearization of A^^ in the Fourier space. The solution u of 
(2) can be expressed as a convolution product of the characteristic function of E 
and of the anisitropic kernel 



However, this kernel (more precisely K^^^^^i) does not satisfy the hypotheses (Ai) 
and {A2) above: K^j^^i is not positive and x -^ J^d |xpi^(^(x) is not in L-^(R). In 
section 2, we etablish some properties of the anisotropic heat kernel X^^. Precisely, 
we prove that the associated Hamiltonian flow is 

which establishes a link between K^f^ and ^-anisotropic mean curvature flow. 

In section 3, we establish the consistency of a Bence-Merriman-Osher scheme 
based on ([5]). We have however not been able to prove the convergence of the 
algorithm to ^-anisotropic mean convergence in the general setting of uniformly 
bounded and continuous functions. The main difficulty in trying to extend the 
argument of [27]. is the thresholding and the lack of monotonicity of our sheme 
that may not preserve the continuity of the front. 

Therefore, in the last section, we present numerical evidence of the conver- 
gence of a modified scheme. In this scheme, the thresholding is obtained via a 
reaction term, in the spirit of phase-field approximation. Computationally, the 
scheme proves very efficient and very fast, even when the anisotropy (j)^ is not 
smooth. 

2 The operator A^^ and properties of the anisotropic 
kernel K^ 

Let (j) — (j){^) denote a strictly convex smooth Finsler metric and let (j)^ denote its 
dual (see [8J). We assume that that (j)^ is a 1-homogenous, symmetric function 
in C^(R^ \ {0}) that satisfies 

A|C|<0^(O<A|C|. (6) 



In particular, it follows that for any ^ G M^ and t G M, 

'rm = \t\r(o 

The associated anisotropic mean curvature is defined as the anisotropic Laplacian 
operator 

A^n = div (<^°(Vn)(/>^(Vn)) , Vn G H^{n) 

A direct computation shows that for any ^ G M*^, 

fA^[cos(27r^.x)] = -47r2(/)°(^)2cos(27r^.x) 
\A^[sin(27r^.x)] = -Air'^cp^'i^f sm{2Tr^.x), 

i.e., that plane waves are eigenfunctions of the anisotropic Laplacian (albeit non- 
linear). We define A4, : H'^{R'^) -^ L'^{R'^) by 

Given an initial condition uq G L^(R^), we study the solution u of, 

\u(0^x) = Uq 

The function u can also be expressed as the convolution product u — K^f^^t * uq^ 
where the anisotropic heat kernel K^j^^t is defined by 



K^^t = ^"' 



-47TH(f)''{0 



We also set K^fy = ^(/),i- In the rest of this section, we establish some properties 
of this operator. 



Proposition 1 (Regularity of K^j^). 

The function K^p • ^ ^ g-^^'^'^^)' is in l^^+^'i 

is a regular function. 



and the distribution D^^^Ka, 



Proof. First, we claim that the Hessian of if^^ is a regular distribution since 



DK^{i) = -8^^ 



%i)e 



and 






of^\\^-4n^-{0\ 



We note that (p^ is discontinuous at ^ = 0. Nevertheless, we next prove that the 
d— 1^^ derivative of D^K(p is a regular distribution, without Dirac mass at ^ = 0. 
Assume that / = D^^'^K(p is an integrable function on R^ for some integer n < d. 
The homogeneity of 0^ shows the existence of a constant Cn such that 

\D-+^K^\ < C„^e-^l«l', for all ^eR^X {0}. 

Since / is smooth away from ^ = 0, the distributional derivative of / is the sum 
of a regular function and of possibly a Dirac mass at ^ = : 

Df = {V/} + c S, 

where c is a constant and V/ denotes the pointwise derivative of /. Let (f G 



V 



^n+2 



and let e > 0. Then 



{Df , (f) = - (/ , div(/9) = - / fAivifda 



f.divipdx — / f.divipdx 



— / Vf.Lpdx— / f.{(p.n)da— / /.divc/^dx, 

jRd\B{0,e) JdB{0,e) JB(0,e) 

Since we assumed that / G L^{R^)^'' , the last integral above tends to 0, as 
e ^ 0. Moreover as n < d, we have 

/ / Lp.nda 

JdB(O.e) 



< Ml- / Cnj^e-^\^\'da 

JdB{0,e) ICr 

< II ^ II L^^C'^ / e~^da < CnW^WL"-^' 

JdB(O.e) 



d—l—n 



JdB(0,e) 



SO that 



lim 



/ / if.ndcr 

JdBiO^e) 



0. 



/dB(0,e) 

It follows that c = 0, which concludes the proof. 
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Proposition 2 (Decay properties of K^fy). 

Let s G [0, 1[. There exists a constant Ccj^o^g^ which only depends on the anisotropy 
(j)^ and on s, such that 

\K4-)\ < 17^^' VxGM^. (7) 

Remark 1. The case s = is easy: According to proposition \l\ the function 
A^~ K(fy{^) is in L^(R^). The continuity of the Fourier transform from L^ to 
L^ shows that 

11(1 + \xf+^)K4L^ < CWK^iO + A^i^^(0||Li(M.), 



and since K(p{^) — e 



i^^(-)i ^ ifSfei' ^-^ 



The proof uses properties of interpolation spaces [TU]. Consider X, Y two 
Banach spaces, and for ix G X + y and t G M^, let 

k{t,u)= inf {||i/o||x + ^||^i||y}- 

For 5 G [0, 1] and p > 1, the interpolation space [X,Y]s^p beetween X and Y is 
defined by 

[X, Y]s,p ^LeX + Y ; t-'K{t, u) G L^ Q 

In particular, given a strictly positive function /i : R^ ^ R, consider the weighted 
space L^ defined by 

L^(M^) = L G L^(R^) ; sup {h{x)u{x)} < oc I . 



One can interpolate between L^ and L^ according to the following lemma. 
Lemma 1. Let h be a strictly positive function R^ -^ R^ and let s g]0, 1[. Then 



Proof. 1) Assume that u G L^(M''). There exists a constant C such that for a.e. 

xeR'^, 

l»WI < ^. (8) 

To estimate k{t,u) — inf„=„Q+„^ \ \\uq\\l°° + iHiiillL"" [, we note that 
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• If t > 1, the choice uq — u and ui — 0, shows that K{t^ u) < \\u\\l<^. 

• If t < 1, we consider the set A = Ix ^R^ ; \u{x)\h{x) < t^~^ k and we 
choose U{) — xa^ u and ui = XA u, so that ||ixi||L°^ ^ t^~^. Moreover, we 
remark that for ah x G A^^ \u(x)\h(x) > t^~^ so that, in view of ([s]) 

\uo{x)\ < Ch{x)-' < C\uo{x)\H'^^-'\ 

and thus k{t, u) < (C ^ l)t\ 

In summary, these estimates show that 

^ ^ ~ \{C + l)t' ift< 1, 

which proves that u G [L^,Lf^]s^oo. 

2) Conversely, we consider u G [L^,L^]s^oo- For ah t > 0, there exists a decom- 
position u = uo^t + ui^t such that 

\uO,t\L^ +t\ui^t\L?^ < Cf. 

15 1 ' ' ' Al 

It fohows that for ah t > 0, we have 

h{xy\u{x)\ < \h{xy \u^^t{x) + ui^t{x)\ < c (h{xyt' + h{xy-H'-^y 

Choosing t — h{x)~^ in the above inequahty shows that for all x G M^, 

h{xy \u{x)\ < 2C, which concludes the proof. D 

We use the following properties of interpolation spaces: 

(Pi) if T is continous from X ^ X and from Y ^ Y^ then T is continous from 

[X,Y]s,pto[X,Y]s,p. 

(P2) if p < p', then [X, Y]s^p C [X, yj^y for any < s < 1 and p>l. 
(P3) [L~(R<^), L^+|,|)(M'^)].,oo = L^+|,|).(K'') for any < s < 1. 

In the foUowing, we consider the case where T is the Fourier transform, X — 
L\R'^), Y = L°°{M.<^), X = W^^\R'^) and Y = i^(?+|^|)( 
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Proof of Proposition^^ We claim that it suffices to show that for any < 5 < 1 

uiO - A^K^(0 e[X,YUi. (9) 

Indeed, the inclusion [X, yj^^i C [X, yj^^oo implies then that u G [X, yj^^oo, so 
that in view of (Pi) and (P3) we obtain 

u e [X,Y]s,oo = [L-(R''),L^+|,|)(M'')],,oo = i^^+|.|)4K'), 
and consequently 

1(1 + \x\')u{x)\ = 1(1 + |x|^+^)K^(x)(l + |x|)"| < C^o^„ for ah x G M^. 
It follows that for some constant C^^o ^ 

We now prove Q. The homogeneity of 4>° shows that for some ci > and 
C2 > 0, and for ^ G M'^ \ {0}, 

l-(OI < ^-""^'^ -d \Vum < |^e-^l«l^ 

which shows that i/ G X = L-^(R^). However, ix may not belong to y = L^{R^). 
We now estimate k{u, t), for t G M^. If t > 1, we set uq = ix, i/i = 0, so that 

k{t,u) < \\u\\x, yt> 1. (10) 

If t < 1, consider the function pt(0 defined by 

'O if \x\<t 

PtiO = I 1 if |x| > 2t 

sin ( I ^^^ j otherwise. 

We choose i/q = (1 — Pt)'^ and ixi = pt'^ and check that 

^ ^ JB{0,2t) JB{0,2t) IQ"^ ^ 

Moreover, 

< /, |VptKOdC+/, |Vi/(0|dC 
2t J B(o,2t)\B(o,t) \t,r JR''\B(o,t) \t,r 

12 



First, we have 



JB(0,2t)\B(0,t) £K-^ ^ - 2i ' ' A - 2 



Second, 



/. 



Rd\i?(0,«) ICI'^ 



-A|CI 



Jb{o,i) 



C 



-A|SP 



l)\B(0,t) ICI'' 

/•I 1 /• 

< C|S<^| / -dr + ClS'^l / 



JM'i\S(0,l 



)W 



-m'- 



d£. 



e"^'' dr 



< CIS*^! |ln(t)| + 



1 a/^ 



\/A 2 



so that 



l^iilly < C 



Tx 1 a/tt 



+ |ln(t)| 



\2 VX 2 
Consequently, this decomposition of u shows that 

k{u,t)<C{l + \ln{t)\)t, Vt<l, 
for some constant C > 0. In summary, 

^ ^ ~ \C(l + |ln(t)|)t if t<l, 
and therefore we obtain 

||t-^fc(t,n)||ii(i/,) = / \k{t,u)t-'\'^dt 

^i(Co + Ci|ln(i)|) 



< 



/ 

Jo 



/o i^ 

which proves that u G [X, l^]s,i as claimed. 



-dt + 



L 



oo |L,||1 



(11) 



(it < +00, 



n 



Corollary 1. For any 5 G [0, 1[ and p G S^^ 



{x(^xK^)^^^(EL'{W'-'). 
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Proposition 3 (Decay of averages of K^f^ on spheres). 
The integral 

i{R) = / K^dn^-\ 

is stricly positive, and decays rapidly as 

where A and A are hounds for (p^ as in ^. 

Proof. Since the measure fi := 5aB(o,i?) has finite mass, its Fourier transform is 
the continuous and bounded function 



jRd Jqi 



—27Tix-^ 



ldB{0,R) 

As fi is radially symmetric, fi can be expressed in the form 

where J is a function R+ -^ R. It follows that 

I{R) = {SdBio,R) , K^) = (i?'-' J(i?ICI) , e-4-'^^(^)' 

= i?^-i / r^ r'-'j(Rr)e-''-'^'^'^'^'drdn'-\ (12) 

We use the particular case when 0^(0 is isotropic, i.e., 0^(0 = ICI to estimate the 

1 -^ 
previous integral. In this case, K(p = 72"4d72 ^ ^ is the heat kernel, and by a direct 

calculation we see that the corresponding integral is I{R) —< 5B{id,R)i^(l) >— 



j^d—l\gd—l\ _R^ I I 

— . L^2 ^ ^ • Comparing this expression to ( 12 ) and using the radial symmetry 
of K^ shows that 

or, after a change of variable, that 

/ rJ(i?r)e-4'^'^''('')''''dr = -— T7^— -.e"J^^W^. (13) 



Returning to a general kernel K^, we deduce from (12) and (13) that 

^ ' (47r)'^/2 7s<*-i <?!'°(6')'^ 

which in view of ([g]) concludes the proof. D 
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Proposition 4 (Positivity on hyperplanes). 
For all p G §^; the integral J ± Kcf)dT-L^~^ is well defined, and satisfies 



[ 



KsdH 



d-l 



2^(P0{p)- 



AA ip-i 



jd-l 



1 



2v^A - Jp± ^ - 2^X 



In particular, we have 



Proof. Let p G S . We already know from Corollary 1 that J ± K^pdM is well 
defined. Consider for // > 0, the approximating functions /^, defined by 

jf^{x) = K^{x)e--\-\y^'\ 



The function /^ belongs to the Schwartz space 5(R ). Moreover, /^ -^ K^fy in 
VI/^~-'^'-'^(R^), and the trace trace theorem [26j shows that one also has 



lim / f^{sp)ds 



K^{sp)ds. 



(14) 



On the other hand, it follows from the Lebesgue dominated convergence theorem 
and from Q that 



lim / fadU"^-^ = / K^dn'^-\ 



(15) 



As /„ G 5 



we infer that 






SO that (14) and (15) yield 



/ K^dU"^-^ = / K^{sp) ds= f e-4^'^'^^(^)' ds 
Jp^ Jr Jr 



-7r(2v^(/)«(p)s 



ds = 



2V^0^(p)' 



which concludes the proof. 



D 
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Proposition 5 (Moments of order 2). 
Let p G §^. Then ^ J ± x xK(fydH^~^ is well defined and satisfies 






2^' 

Proof. Corollary nl states that the integral / i \x\^K(i)dT-L^~^ is well defined. Re- 
calling the sequence /^ used in the previous proposition, we observe that i?^//i -^ 
D^K(f) in VI/^~^'^(R^), so that the trace theorem implies 



lim [ D^f^(sp)ds = f D^K^(sp)ds. 

From proposition [2] and the Lebesque dominated convergence, we obtain 

lim / x®x fJx) dU"^-^ -^ / x®x KJx) dW 
^i^ooJp± ^ Jp± 



jd-i 



(16) 



(17) 



Moreover, we have 



jd-i 



^x®xf^{x)dn = \5p^ , x®xf^J ^ -^{5p, D f^ 



—^jDUsp)ds, 



SO that in view of (16) 



/ x®xKJx)dH^-^ ^-—^ I D^KJsp)ds. 

Jp^ 47r^ Jr 

We next estimate the above right-hand side by a direct calculation: 
-l^f D^K^{sp)ds = [24>''{p)(t>h{p) + 2ct>l{p)®(t>l{p)\l e-^^'^'^°^P^'ds 

^'* J R J R 

- '2mp)®mp)\ [87Th^%pfe-'^^''''^^^P^'ds. 

^ ^ ^ ^ Jr 

Further, we see by integration by parts that 

/ Snh^ripfe-^^"'"^"'^^^" ds = [ |47r22s</.°(p)2e-4-'*'*''(p)' | {s} ds 

Jr 2^(po{p) 

16 



and we conclude that 
1 

2 



X xK^{x)d'H 



d-l 



<Ph{p) 



2v^' 
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Corollary 2 (The operator F{X,p)). 
Given X eR'^'"^ and p e S<^, let 



F{X,p) = ([ K^{x)dn'^-A Q / <Xx,x> K^{x)dH'^-A . (18) 



This operator is elliptic and satisfies 

FiX,p) = ci>°{p)4>l^{p):X. 



(19) 



Proof. Equation (19) is a direct consequence of propositions [4] and [5J while the 
elhpticity of F fohows from the convexity of 0^. D 

Remark 2. In the next section, we introduce an algorithm for motion by anisotropic 

mean curvature, and show its consistency with an evolution equation of the form 

\/u [ I [ I 

ut = —FiD'^u, ) where F is defined by (18). The expression (19) shows that 

\Vu\ 



this operator is precisely the one corresponding to motion by anisotropic mean 
curvature (see f^). 



Proposition 6 (Positivity of order moment s). Let V be a subspace of R.^ of 
dimension 1 < m < d, and let < s < 2. Then 



b 



\x\'Ksdn'^ > 0. 



Proof. We first consider the case m — d and V — W^. we consider the finite part 
Pf ( I }d+s ) ^s a temperate distribution, defined for Lp G 5(R^) by 



pf 



\d+s 



(f) = hm 



lim < / 



(f{x) - (p{0) 



Rd\B(0,e) \X 



d+s 



dx 



This function happens to be the Fourier transform of the distribution |x|^. More 
precisely, 



J'Wxi 



Cs^Pf 



1 



12^1 



d+s 



with C., = 2-+V/^ ^^^^ + y ,(20) 

s,d r(-s/2) ' ^ ^ 
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(see for instance ^25j, F denotes the Gamma function). We can thus write 

lj-\'K,dx = (1x1% K,) = (C,.P/ (^^) , e-^-^^^^^^^ (21) 

= ^^'^ H / . 1^ .1^+. > 0^ (22) 

a stricly positive quantity, in view of the sign of Cs^d- 

Suppose now that m < d and consider the subspace V = Vectjei, . . . , e^}. 
We write x = (x',x'0, C = (C', ^^''''^''0. with x',f G y. A straightforward 
computation shows that 

/ \x I Kcfydri = (|x | , K^f^yx ^^))x>'(Mjn)^v(w^) 

c^^,-^/(^^^)'Mn) 

where the function h : W^ ^ R is defined by 



The next lemma states that h is C^ and maximal at ^^ = 0, which in view of (20) 



and of the sign of Cg^m concludes the proof. D 

Lemma 2. The function h : R^ -^ R, defined by 

is C^ , with fast decay as \^^\ ^ oo, and is maximal at ^' — {). 
Proof, recalling ([g]), we first remark that 

so that the functions ^' -^ e~^^ ^""^^ '^^^^ and their derivatives are uniformly 
bounded in L^{R.^~^). The C^ regularity of h is thus a consequence of the 
Lebesgue theorem. The above estimate also shows that 
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To determine the maximal value of h, we consider the sets A^/^j, defined for 
alU' GM*" and t g]0, l[by 

A^,^t = {f e M-^— ; e-4-''^''(«''«"))' > t} 



M,^.' 



% 


Br. 


y' 


k%,% 


Aol yV 




V 


I, 



Figure 1: 



Fix ^0 ^ IK^- The set A^ ^ can be defined as the intersection of the hyperplane 
|C ^ M"^ ; f = Co} with the Frank shape 



B^o^, = |C G R^ 0^(0 < ^7-/n(t)| 



The set B(^o^i is convex since 0^ is convex. Moreover, from the symmetry of 0^, 
('/'°(0 = </'°(-6),wehave 



1^ 



€i>tl 



|A 



■So.*l- 



Next, let 

We remark that the convexity of (\f implies that A^ ^t C Ao,t- Indeed, let 



</'°((0,O) = </'''Q((e^,el') + (-Co,C2)) 



< 2 ('^° ((^o,er)) + r ((-e^,e^'))) < ^\/-^"(^)' 
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so that e-^^ ^"((o^^'O) > t, i.e. f G Ao^f Invoking the Brunn-Minkowski in- 
equahty, we obtain 

lAej'^^"-""^ = l\A^'^,t + A_^,//^^-"^^ (23) 

> i (|A^,,,|V('^— ) + |A_^,,//('^— )) > \A^,//^^-"^\ (24) 

and finally that, 

l^o,t| > \A^^J > \A^^J- 

As this equality holds for any ^q ^ ^^^ it follows that /i is maximal at ^^ = 0. D 

3 The Bence-Merriman-Osher-like algorithm 

Barles and Souganidis [H] have studied the convergence of a general approximation 
scheme to viscosity solutions of nonlinear second-order parabolic PDE's of the 
type 

ut + F{D^u,Du) = 0. (25) 

The main assumption on the function F is its ellipticity, i.e., F satisfies 

VpGM^\{0},VX,y GMf ^, X<Y ^ F{X,p) <F{Y,p). (26) 

Barles and Souganidis study a family of operators Gh '■ BUC{R^) -^ BUC{R^) 
for /i > 0, which satisfy, for ah u,v e BUC{R^) 

• Continuity 

VcgM, Gh{u + c) = GhU, (27) 

• Monotonicity 

u<v ^ GhU < GhV + o{h) (28) 

(see remark 2.1 in [6J) 

• Consistency 
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For all T > and for all partitions P = {O = to < ... < tn = T} of [0,T], one 
can then define a sequence of fonctions np : M*^ x [0, T] — )■ R by 

fG,_,,(np(.,t.)) if te(t.,t.+i], 

[g if t = 0, 

If additionnally the following condition holds, 
• Stability 

there exists u G C([0, oc], [0, oc]), independent of P and depending 
on g only through the modulus of continuity of g^ 
such that Cl;(0) = and for all t G [0,t], 

J\up{.,t) -^||l- < ^(^), 

then the following theorem holds [5] : 

Theorem 1. Assume that Gh : SC/C7(M^) ^ St/C(R^) satisfies ^, ^, ^, 



and ^ for allT > 0, g e BUCiR"^) and all partitions P of [0,T]. Then, up 



defined in (30) converges uniformly in R x [0,T] to the viscosity solution of (25). 



This result was used by H. Ishii, G. Pires and P.E. Souganidis in [27j to study 
anisotropic mean curvature flow. These authors introduce a kernel /, which 
satisfies: 

(Hi) f{x) > 0, f{-x) = f{x) for ah x G R^, and J^d f{x)dx = 1 

{H2) jp±{l^\x\^)\f{x)\dn'^-^ <oo for ah p ^ S^ 

J the functions p -^ J^± f{x)dH^~^ p -^ /^± XiXjf{x)dH^~^, 
[^ ^ hj ^ d^ are continuous on S^ 

{H4) J^d |xp|/(x)|dx < 00 

(i^s) For all collections {R{p)}o<p<i C R such that R{p) -^ oc and pR{p)^ -^ 
as p ^ 0, and for all functions g : R^"-*^ ^ R of the form g{^) = 
a + iA^, with a G R and A G S^-\ 



lim sup sup 



Jb(0,R(p)) JRd-i 



0, 



where 0(n) denotes the group of d x d orthogonal matrices, and where 

/[/ : R-^ ^ R is defined for all U £ 0{d) by /[/(x) = f{U*x). 
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Theorem [T] has been applied to schemes for anisotropic mean curvature motion 
(see theorem 3.3 in [27]) with Gh defined by 

G/,*(x) = sup{A e R ; ShlM,>xix) > Oh} (32) 

= mf{XeR; Shl^>xix)<eh} (33) 
where 

Sh9{x) = h-'"^f{./Vh) * g{x) = h-^/'' [ fiy/Vh)g{x - y) dy, Oh ^ ^ + cVh, 

and where F{X^p) is given by 

F{X,p) = - l^j ^ f{x)dn''-\x)^ Q j ^ {Xx , x) f{x)dn''-\x) + c|p|) , 

(the last term in this integral models a forcing term). 

In this section, we follow the proof in [27j to show a consistency result in our 
case when f is a non positive kernel and does not have moments of order two ( 
ie. X -^ |xp/(x) ^ L^iW^)). We introduce two operators G^ and G^ defined by 

G+^{x) = sup{A G M ; Sht^>x{x) > Oh} (34) 

Gl^{x) = inf {A G M ; Sht^>x{x) < Oh] (35) 

which are not necessarly equal as our kernel is not being nonnegative. 

3.1 A consistency result in the case where / = K^j) 

To adapt these results to our context we modify the assumptions (iifi), (^4) and 
(i^s) as follows 

{H[) Jp^ f{x)dn^-^ > for all p G S^, f{-x) = f{x) and j^d f{x)dx = 1, 

{H'^) j^a \x\^-^\f{x)\dx < 00 for < /x < 2, 

{H'^) Assume that ji g]0, 1/2]. Then for all collections {i?(p)}o<p<i C M such that 
R{p) -^ 00 and pR^pf'^ ^ as p ^ 0, and for ah functions g : R^"^ -^ R 
of the form g[i) = a + (A^ , ^ with a G M and A G S^-\ 



lim sup sup 

P^^ UeO{d) 0<r<p 



b(0Mp)) 



lim sup sup 

P^^UeO(d)0<r<p 



jRd-l 

Jb(o,r(p)) JR"*-! 



0, 



In this last statement, B{0,R{p)) denotes the (n — 1) -dimensional ball, cen- 
tered at and of radius R{p). 
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3.1.1 K^ satisfies {H2,H^) and {H[,H'^,H'^) 

We remark that k^{i) = K^{-i) and J'(K^)(0) = 1, so that 

KA—x) — K(f){x) for all x G M , and / KA){x)dx — 1. 

Moreover, proposition Q shows that 

j^^K,{x)dn'-^>j^;^,>^ for all p e S\ 

SO that {H'l) is satisfied. Propositions Q and ([5]) also imply that K^j) satisfies 
{H2), i.e., 

/ {l + \x\^)\K4,{x)\dn'^-^ < 00 for all peS'^. (36) 



Concerning (H^), we note that 




I [ X ® xK^{x)dH'^-^ 


= 2^«<'''' 


and that 






1 


2v^</.°(p) ■ 



Since (j)° is smooth on W^ \ {0} and positive (in particular (j)° > \ on^'^ ) we see 
that the functions 



P 



f K4,{x)dn'^~^ P^ [ XiXjK^{x)dn'^~'^, l<i,j<d, 



are continuous on S^. 

We next prove that if < /x < 2, then 

/ \x\'^~^\f{x)\dx < oc. 

JRd 

Indeed, proposition [2] with 5 = 1 — fi/2 shows that 

\x\^~^\f(x)\dx < f ^r,s\x\ , r c 



roo 1 

< CIS'^I / T^TT^dr < oo, 



(1 + ;.l+A»/2) 
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for some generic constant C. It remains to prove {H'^)\ Let < /i < 1/2 and let 
R : R+ — ^ R+ such that, as p ^ 0, R{p) -^ oc and pR{p)'^-^ -^ 0. Setting 
fu{x) = K(fy{U'^x)^ we consider 



/ fu{x,rg{x))g{x)dx^ / . 2-^\ fu{x,rg{x))g{x)dx 

Jb{o,r{p)) Jb I oMp) ^~ ) 



Let /ip, /i : 



+ / / 2-^.\ fu{x,rg{x))g{x)di. 

Jb{oMp))\b(oMp)^) 

■ R denote the functions 

hpr(x) = fu(x,rg(x))g(x)x 2-/^ 



(37) 
(38) 



h{x) = fu{x,0)g{x). 
When r < p^ hp^r{x) converge to h{x) pointwise as p ^ 0, and 



\hp,r{x)\ < 



C 



1 + 1x1^-1+^ 



gl 



1/TD)C^-1\ 



for some constant C independent of p, r and t/. Invoking the Lebesque dominated 
convergence theorem, we conclude that 



lim 

p— ^0, r<p 



/ hpr{x) dx ^ / h{x)dx, 



uniformly with respect to U and r. The second term in (37) converges to 
uniformly with respect to U and r as p ^ 0, since 



/. 



b{o,r{p))\b(oMp)^) 



\fu{x,rg{x))g{x)\dx 



< C 



L 



b(o,r(p))\b(o,r(p)^^ 1 + 1^ 



d-l+s 



dx < CIS 



d- 



1 f^(p) !_ 

/ 2-/X :; . r~ 
jR{p)^r- 1 + |r 



1+5 



dr 



< CIS 



d-i\ 



R{p) 



-(2-/x)s 



R{p)- 



for some generic constant C. We conclude that 



lim sup sup 

P^^ UeO{d) 0<r<p 



Jb\ 



fu{x,rg{x))g{x)dx- / fu{x,0)g{x) dx 



ib(0Mp)) 
The second statement in (iJg) is established similarly. 



0. 
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3.1.2 The consistency proof 



Proposition 7. Let (f G C^(R^). For all z e W^ and e > 0, there exists 6 > 
such that for all x G B{z^ 5) and h G (0, 5], if V(j){x) ^ Q we have 

Gj^(f{x) < if{x)^{-F{D^if{z),Dif{z)) + e)h 
and G^(f(x) > (f(x) + {-F{D^(f, D(f(z)) - e)h. 

Proof We closely follow the argument in [27j . 

1. We only prove the first inequality The other one is obtained similarly 

2. Without loss of generality we can assume that z = 0. Let us fix a G M, 
such that 

a> -F{D^(p{0),D(p{0)). 

The inequality is proved if we can exhibit a 5 > such that, for all x G 5(0, S) 

and h G (0, 5], 

Sh\>ip{x)^ah{x) < ^h- 

3. Fix 6i > 0, such that D(f 7^ on B(0,6i) and choose a continuous family 
{U{x)}x^B{o,Si) ^ 0{d)^ such that for all x G 5(0, 5i), 



U{x, 



y\Dy,ix) 



where ed denotes the unit vector with components (0, 0, ..., 0, 1) G M^. Note that 
if X G 5(0, 5i), then 

4. Choosing 5 smaller if necessary, (5{) implies the inequality 

a > -F{D^(f, Dip) in 5(0, 5i), 
or in other words, 

< -c\D^ix)\, (39) 

where P denotes the d x (d — 1) matrix with components Pij = 6ij. 
5. We next fix e > 0, and S2 G (0, 5i[, such that for all x G 5(0, ^2), 

^^_^ (P*f/(0)(I)V(0) + 3e2/)f/(0)*PC, e)/t/(x)(C,0)d^ 

-{a-e^)f fu^,){^,0)dC< -{C + e)\D^m. (40) 
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6. The Taylor theorem yields a 7 > such that for all h > 0, y e M. , and 
X e B{0, S2), if \^\y\ < 7, then 

(p{x - VhU{x)*y) < (p{x) - Vh {Dip{x) , U{xyy) 

+^{uix)iD\ix) + e^I)U{xry,y 

< <fix) - Vh\D<f{x)\yd + Chyj 

+^ {p*Uix)iDM^) + 2e^I)UixrPy' , y') , 



and 



<f{x - VhU{x)*y) > ^(x) - Vh {D(f{x) , U{x)*y) 

+^(u{x){DMx)-eH)U{xry,y 
> <f{x) - Vh\Dif{x) \yd - Chyj 



+- {P*U{x){D^^{x) - 2e^I)U{xyPy' , y') , 

where we write y — {y\ yd) G W^~^ x R, and where C is a positive constant. 
7. Reducing 7 and 82 if necessary, the previous inequahties imply that for 

y G 5(0, 7/ V/^) and x G 5(0,52), 

• if ip{x — vhU{xyy) > (p{x) + a/i, then 
Vh 



Vd < 
< 

. if 

Vd < 
then 

We define 



Vh 



Vh 



{-a + e' + ^ (P*f/(0)(DV(0) + 3e2/)f/(0)*Py' , y' 



\D^iO)\ V "''"''+ 2 (^*f^(0)(^V(0) - 3e^I)Ui0rPy' , y' 



(y?(a: — VhU{x)*y) > (p{x) + ah. 



a' ^ {a - e'^)\Dip{0)\-^ 

a, = {a + e^)\D^m-' 

A' = \D^{0)\-^P*U{0) (DV(0) + 3e2/) t/(0)*P 

A, = |D</?(0)|-ip*f/(0) (r>V(0) - 3e2/) C/(0)*P, 
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and for y' G M'^"^ 

9'{y') = [-a' + \ {A'y' , y')) 9e{y') = (-a. + ^ {A,y' , y')) . 

We also set 

Vh,x ^lyeR"^; (f{x - VhU{xyy) > ^{x) + a/i} , 
and 

'Kh,. = {y e K'^ ; yd < Vhg,{y') 
E-h,.-{y^^'^-^ yd<Vhg'{y')\. 

We check that for all x £ B{0, S2), 



[e-^^^ n 5(0, 7/v^)) c (y,,,, n 5(o, 7/ v^)^ 



8. The assumption {H^) yields the existence of a decreasing function uj G 
C([0, 00), [0, 00)) such that uj{R) ^ as i? ^ oo, and 

/ |/(y)||y|2-'^dy < a;(i?)2, for all i? > 0. 

JB(0,fi)<= 

For each < t < 1, we define the family of sets R{t) G (0, oo) by 

uj{R{t)) = tRitf-^, (41) 

which satify (-^5). We then choose r £ (0, 1) such that 

R{t) < l/t, for all t e (0, r] (42) 

9. Let 

p=Vh, T{p) = S„_i(0, R{p)) X M C M"*. 
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< 



For all h g]0, t^) and for all x G 5(0, ^2), we estimate 

/ fu(x)(y)dy = / fu(x)(y)'^^>^(x)+ah(x - VhU%x)y)dy 

/ fu{x){y)dy+ / \fu{x){y)\dy 

< / _, fu{x){x)dx+ \fu(x){y)\dy 

+ L ^ . \fu{x){y)\dy 

< /^ fu(x){y)dy+ \fu{x){y)\dy 
+3 / \fuix)iy)\dy 

Jb(0,R(p)Y 



'B(0,R(P)) 

10. For the last integral above, we have 



/ \fu(x)\iy)dy < -7^2^ I \y?-''\fuix)\{y)dy < oj{R{p))p, 



and moreover, since K^ is symmetric, 
1 



o / fu{x){y)dy < / fu{x)\{y)dy + uj{Rip))p. 

We note that 

/ _, fu{x){y)dy^ / fu(x){y)dy 

JT(p)nE+^^^ JT(p)n{ya<pgHy')} 

f r fpg^iy') 

= fu{x){y)dy+ d^ fu{x)i^,r)dr 

= / fu(x){y)dy+ dr fuix){C,rg{^))g''{^)d^. 

JT(p)n{ya<0} Jo J B„-i(0,R{p)) 

It follows from (-^5) that as p — )■ 0, 

-{ ^ fu{x){y)dy- fuix){y)dy\ -^ / /c/(x)(^,0)/(Od^, 

uniformly with respect to x. Possibly reducing r we may assume that for x G 

5(0,52), 

-If , fu(x){y)dy- f fuix){y)dy\ < I /c/(x)(C,o)/(CK + e' 
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Using same argument, we also conclude that 

/ ,^ Afu{x){y)\dy = \i \fu(x){y)\dy\ 

-(/ \fuio:){y)\dy] 

[jTip)n{0<yd<P9.{y')} J 

< Pf^ \fuix)\(mi9H0-9emd^ + pe 

< pe^ (l + j^^_^ (2 + 3|C|2) |/t;(,) |(C, mi 

< Cope\ 

where 

Co= sup |l+/ (2 + m^)\fui.)\{i,0)dA. 



-c- e, 



11. Finally, noting that from (39), 

/ fui.)ii,0)9'i0di < 
we get 

/ f{x)\><p{x)+ah{x-'/hz)dz < -+/ fu{x){C,0)g{C)d^ 

+p (e^ + 4uj{R{p)) + Coe 

< -+p(-c-e + e2 + 4a;(i?(p)) + Coe' 

< ^h, 

for e sufficiently small. D 

Even if the function (j) is regular, G^(^ and G/^^ need not be equal and 
continuous. However, it is easy to check that if (/:? = 1^ is a characteristic function 
then G^Iq = Gj^^^- The next proposition shows that if (f is smooth, G^Lp{x) = 
G^{x)Lp^o{h)^ so that one could conceivably build a Bence Merriman Osher type 
scheme using either G^ or G^. 



Proposition 8. Let Lp G C'^iW^). Let x eW^ such as V(f{x) ^ 0, then 

G^(/p(x) = G+(/p(x) + o(/i). 
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proof. Let x G M^ such as V(p{x) ^ and for all /i > let 

Introduce also ^/^(A) : R ^ R defined by 

gh{X) = ShX^>x{x) = / K^,h{y)x^>x{x - y)dy. 

This function may not be continuous. We claim that its jumps are bounded by 
o{y/h). Indeed, for all A G R, one can express ^/^(A) as 

9hW = / K^Mxw>x}{x-y)dy+ / K^Mxw>x}{x - y)dy 



where a is chosen sufficiently small so that \V(f{y)\ > for all y G B{x,a). Let 
< /i < 1, let 

^{R)- f \y\^-^\K^{y)\dy, 

and let R{t) be defined by the equality oj{R{t)) = tR{tf-^. Note that (i7^) 
implies that \fhR{hY~^ ^ as /i ^ 0, so that \fhR{hY~^/'^ < a ioi h sufficiently 
small, it follows that 

\Rh{x)\ < [ ^ ^ \K^My)\dy- 

Moreover, changing variables, we see that 
/ ^ \K4,,h{y)\dy < I \K4y)\dy 

^ p^M(2 u)y2 I , \y?~^\KAy)\dy 

< u;(i?(/i)) 



i?(/i)(2-/^)V2- 

Since < (2 - /i)/2 < 1, it follows that 

Further, the fact that |V(/^(y)| > on B(x^a) show that ^/^ is continuous in 
A, which proves the claim. 
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Recall that 

(GJ;^{x) = inf {5 G M ; ShX^>s{x) < 9h} 
\g^(P{x) = sup {5 G M ; ShX^>s{x) > 9h} , 

it follows from the claim above that 

^hX^>G-^(x)(^) = ^h + o{Vh), and 5/^X^>G+^(x)(^) ^^h + o{Vh), 
and consequently 

One can use the same argument as in the consistency proof, (in particular see 
point 7) to show that asymptotically, the above integral behaves like 

where, p = \^'ll(% - In conclusion, as J ± Krj,{x)d'H'^~^{x) > 0, we deduce that 

|V^(x)| 

which proves the proposition. D 

3.2 Discution 

Our consistency result sheds light on the relationship between the kernel K^ and 
the Hamilton Jacobi equation Q. Proving convergence of a Bence Merriman Os- 
her type algorithm in our context seems to be very difficult (if true at all). The 
argument of [27] does not apply here. The main difficulty is that G^ (f may not 
be continuous, even if (p is regular. Further, we can only show monotonicity of 
the operators G^ up to o(h) for smooth functions whose gradients do not vanish. 
The source of these difficulties is really the thresholding in the definition of G^ . 

Thus, rather than advocating for a BMO algorithm, we have considered in the 
next section a numerical scheme where instead of this thresholding, we modify 
the convolution product K^jy^^^Lp using a reaction operator, in the spirit of a phase 
field algorithm. More precisely, given a small parameter e > 0, we may define 
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where T^^^ is defined as follows: Given A G M, r^,e(A) = ip{X) where ip is the 
solution of the ODE 

U(0) =A, 

and W a double well potential with wells located at '0 = and '0 = 1. Note that 
if (/? = 1^ is a characteristic function, then 

4 Numerical simulations 

In the previous section, we proved a consistency result for a Bence Merriman 
Osher-type algorithm. Here we numerically investigate the convergence proper- 
ties of a related scheme, based on a phase- field discretization. We explained above 
why we did not directly implement a BMA algorithm. In the next paragraph, we 
describe the phase- field algorithm based on the operator A^^. 

4.1 The A(^-phase field model and its discretisation 

As an approximation to the anisotropic Allen-Cahn equation Q, we consider the 
following phase-field model 

We also report tests, where we estimate the L^-error on anisotropic Wulff sets 
(the sets which minimize the anisotropic perimeter under a volume constraint). 
To impose volume conservation, we consider a conserved phase-field model, of 
the form 

The parameter 

J^dW\u{x,t))dx 



m 



^/2W{u{x,t))dx' 



can be seen as a Langange multiplier, which preserves the mass of u. See [llj 
where schemes of this form have been studied for isotropic mean curvature with 
a volume constraint. 
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We now describe the numerical method we use for solving the PDE's (43) and 



(44). Several studies of classical numerical schemes for the Allen-Cahn equation 
have already been conducted in the past: see for instance, (THl EH [T31 [El E21 
[SUES]- Here, the computational domain is the fixed box Q = [—1/2,1/2]^ C 
M^, d = 2,3. The initial datum is uq = ^( dist(x,a^o >^^ ^^^^^ ^^ .^ ^ gj^ooth 
bounded set strictly contained Q. We assume that during the evolution, the 
set ^e{t) '-— {ue(x^t) = 1/2} remains strictly inside Q, so that we may impose 
periodic boundary conditions on dQ. 

Our strategy consists in representing ix as a Fourier series in Q, and in using a 
splitting method. First, one applies the diffusion operator, which given the form 
of A(^, merely amounts to a multiplication in the Fourier space. The interesting 
feature of our approach is that this step is fast and very accurate. Next, the 
reaction term is applied. 

More precisely, u^{x^tn) at time t^ = to + n5t is approximated by 

maxi<i<d \pi\<P 

In the diffusion step, we set 

nf (x, t„ + 1/2) = Y. ^..,p(tn)e-4'^'^*^°(^)'e2-P-. 

maxi<i<d \pi\<P 

We then integrate the reaction terms 

U^{x,tn + l) = U^ {X,tn + I) -5te^Wl^{u^{x,tn^ 1/2)). 

In practice, the first step is performed via a fast Fourier transform, with a com- 
putational cost 0(P^ln(P)). 



The corresponding numerical scheme turns out to be stable when solving (43), 

under the condition 5t < Me^, where M = sup^^[o,i] i ^''(^) f • Numeri- 
cally, we observed that this condition is also sufficient for the conserved potential 
in (44). In the simulations, we used W{s) = ^5^(1 — 5)^. 

The isotropic version of our splitting scheme has been studied in ^llj. It is 
shown there that this scheme converges with the same rate as phase-field ap- 
proximations based on a spatial discretization by finite differences or by finite 
elements. Its advantages are greater precision, and unconditionnal stability. 
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Figure 2: Wulff Set (blue) and Frank diagram (red) for the anisotropic densities 
(01,0?), ((/)2,(/)^)and((/)3,(/)g) 



4.2 Test of convergence in dimension 2 

We consider following anisotropic densities 






(161^ 



(161^ + 16 



001 



ki + ^6r°^ + iki-^6 



VSc 



11,001 \ i>ooi 



See figure ( |4.2[ ) for a representation of their Wulff sets B(p. and Frank diagrams 

1. Evolution from a Wulff set. 
We consider the equation 

(dtu = K^u - ^w'{u) 

[^/(O, x) — q (dist{x^ Q.o)/e^) , 
where the initial set f^o is a Wulff set of radius i?o = 0.25 

J^o = |x gM^ ; (/)(x) < i?o}. 

It is well known that the set Vt{t) obtained from f^o through evolution by anisotropic 
mean curvature is a Wulff set with radius R{t) = \/ ^o ~ ^t, which decreases to 

a point at the extinction time text — ^- I^ these simulations, the number of 
Fourier modes is P = 2^, and the time step and phase-field parameter are cho- 



sen to be 5t — \jP^ and e = 1/P. On figure (4.2) the interface Vl{t) is plotted 
at different times. We observe a good agreement between the theoretical and 
computed curves, in spite of the smoothening of the corners of the latter. 
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Figure 3: Q{t) at different times for the anisotropic densities 0i,02,03 



2. Convergence to the Wulff set 



This smoothening of corners actually depends on the thickness e of the diffuse 
interface, as evidenced in the next series of tests, of evolution by anisotropic mean 

The initial set rio is a 

i?o}. 



curvature under a volume constraint according to (44) 

r 



The 



circle centered at 0, of the same volume as fi* = < x G M ; (j){x) < 
evolution Q{t) from Qq is expected to converge to the Wulff set fi*. 

Figures |4.2| -a,b represent the final sets fi* obtained from the resolution of 
anisotropic Allen-Cahn equation, with respective anisotropic densities ^i and 
02, and for differents value of e. We observe that the smaller 6, the better the 

the L^ 



approximation of the Wulff set. In figure 



4.2 



error 



is plotted in a logarithmic scale. This graph indicates that this error is of order 
e. 



4.3 Some 3D simulations 

As final illustrations, we consider the anisotropic densities 




161 + 161 + 161 

The corresponding Wulff sets and Frank diagrams are plotted in figure ([s]). 

We report in figure ^ (respectively in figure ^) the evolution by (j)^ (resp. 
^5) anisotropic mean curvature from an initial torus. The number of Fourier 
modes is P = 2^, the time step and diffuse interface thickness are 6t = 1/P^ and 
6 = 1/P. 
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Figure 4: From left to right : ^{t) at different times with anisotropy ^f , ft{t) 
at different times with anisotropy (/)2, error estimate e -^ W'^b^ — ^B^ r IIli(r^) i^ 
logarithmic scale ( ^f in red and 02 in blue) 



4> 






Figure 5: Franck diagramm and Wulff set : S^^o, B^j^^ , S^^o, B^j^^ 








Figure 6: (/)4(^)-evolution from an initial torus, at different times 
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Figure 7: (/)5(^)-evolution from an initial torus, at different times 
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